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Abstract. We investigate the emissivity properties of the main near-IR transitions of the Fe + ion in the conditions 
prevailing in the inner regions of jets from young stars, based on a simplified 16-level atom model. We present 
new diagnostic diagrams involving prominent near-IR line ratios that allow us to constrain the electronic density, 
temperature, and Fe gas phase abundance ratio, independently of the heating process. Comparison with recent 
near-IR observations of a sample of HH objects indicates gas phase Fe abundances ranging from 15-50 % up to 
100 % of the solar value (depending on the assumed temperature and on the HH object) in agreement with the 
moderate depletions previously derived from optical line ratios or shock models. Hence, it appears that Fe-bearing 
dust is efficiently destroyed in stellar jets. We then use our Fe + emissivity model to predict near-IR [Fe n] emission 
maps for self-similar, cold MHD disk wind models. We show that near-IR [Fe n] lines are stronger than [S n] A6731 
and [O i] A6300 in the cool regions (T < 7000 K) near the wind base, and that observations in [Fe n] with AMBER 
on the VLTI could severely constrain the MHD solution and the inner launch radius of the jet. We also compare 
theoretical predictions with recent observations in the near-IR [Fe n] lines of the L1551-IRS5 and DG Tau jets. 
The cold disk wind model reproduces quite well the two velocity components observed at —100 and —300 kms -1 , 
although the high velocity component appears overestimated by a factor of 1.5 in the DG Tau jet. However, 
the model predicts too little emission at intermediate-velocity and insufficient densities. Similar problems were 
encountered in previous model comparisons in the optical range with jets from T Tauri stars. Denser disk winds 
with stronger heating at the jet base, which have been invoked for optical jets, also appear needed in younger, 
embedded Class I jet sources. 
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1. Introduction 

Constraining the physical mechanism by which mass is 
ejected from young accreting stars and collimated into 
jets remains a fundamental open problem in star forma- 
tion theory. The fact that a correlation between accretion 
and ejection signatures is seen from the embedded Class 
and Class I phases (ages a few 10 4 -10 5 yr) down to the 
optically revealed T Tauri phase (at ages a few 10 6 yr), 
with a similar ejection to accretion ratio, suggests that 
the same mass-loss mechanism may be at play through- 
out the different stages of star formation (see e.g. Cabrit 
2002). Studies of the inner jet structure at each phase of 
evolution are required to confirm the validity of this sce- 
nario. Attempts to compare observations and MHD wind 
predictions have first concentrated on the small scale jets 
associated with optically revealed T Tauri stars (Cabrit et 
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al.Hnm Garcia et al. l2001al2001bl Lavalley-Fouquet et al. 
127)07)1 Bacciotti et al. l27)07)l Shang et al. 1199812002) . These 
studies are conducted in the strong optical forbidden lines 
of [O i] A6300, [S n] A6731 and [N n] A6584. Thanks to 
recent advances in near-IR instrumentation, similar work 
can now be extended to the inner regions of jets associated 
with the younger more extincted Class I sources. 

The imaging study of embedded jet sources with 
NICMOS/HST made by Reipurth et al. <|27jrjTj|l clearly 
demonstrated the importance of [Fe n] near-IR lines as a 
tracer of the innermost jet regions, less affected by extinc- 
tion than their optical counterparts. Furthermore, Nisini 
et al. H2002fl recently performed a near-IR spectroscopic 
survey (from 0.95 to 2.5 (im) of a sample of HH jets, pro- 
viding a first database to investigate the excitation condi- 
tions of the Fe + ion in young stellar jets. Detailed mod- 
eling is now needed to see how such observations can be 
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Fig. 1. Energy level diagram with the prominent opti- 
cal and near-infrared [Fe n] emission lines (wavelengths 
in /jm) originating from the 16 first low-lying levels. We 
also indicate on the left the range of corresponding upper 
level energies (E/k expressed in K). 
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Fig. 2. Theoretical line ratios as a function of n e for differ- 
ent temperature values T c (3000, 10 000, 20 000 K). The 
curves are deduced from our 16-level model and symbols 
from a 142-level computation (Pradhan, private communi- 
cation) . T e increases from bottom to top in panels (b) , (c) 
and (d), from top to bottom in panel (a). 



used to constrain the jet driving mechanism in embedded 
sources. 

We investigate in this article the emissivity properties 
of the main near-IR lines of the Fe + ion, based on a sim- 
plified 16-level model, in the conditions prevailing in the 
inner regions of jets from young stars. Extending the work 
of Nisini et al. l|2TiU2"|l and Beck-Winchatz et al. GUUU, 
we first derive in Sect. 2 diagnostic diagrams that allow to 
constrain, independently of the heating process and ejec- 
tion model, the electronic temperature and density as well 
as the Fe + /S + gas phase abundance ratio — a critical pa- 
rameter to derive mass loss rates from [Fe n] line fluxes and 
to estimate the amount of depletion onto dust grains. We 
then turn in Sect. 3 to detailed predictions for a specific 
class of ejection model, the self-similar cold disk wind so- 
lutions developed by Ferreira l|1997fl . for which consistent 
thermal solutions have been recently computed by Garcia 
et al. I|2001afl . Our study extends the work of Garcia et 
al. H2001b|l . which was restricted to the optical domain. 
We compare these predictions with recent observations in 
near-IR [Fe n] lines of young stellar jets from embedded 
sources, in particular the velocity resolved spectroscopy 
of the inner regions of the L1551-IRS5 (Pyo et al. IMJ^ 
and the DG Tau jets (Pyo et al. l2003|l . We summarize our 
findings in Sect. 4. 



2. Near-IR [Fe n] diagnostic diagrams 

2.1. The 16-level Fe + model 

To predict near-IR [Fe Ii] line emissivities, we build a 
NLTE model including the 16 first fine structure lev- 
els, which comprise the 4 terms (3d 6 4s) 6 D, (3d 7 ) 4 F, 
(3d 6 4s) 4 D, (3d 7 ) 4 P (Fig. Q}. Levels in the 4 F multiplet 
lie ~ 1.7 eV above the ground state so that our sim- 
plified model will be valid for excitation temperatures 
below 2 x 10 4 K, a condition typically achieved in stel- 
lar jets (see below). All radiative transitions among these 
four terms are of quadrupole or magnetic dipole type. We 
therefore assume optically thin conditions and compute 
level populations under statistical equilibrium including 
electron collisional excitation and spontaneous radiative 
emission processes. We take transition probabilities from 
Nussbaumer & Storey 1|1988[1 . level energies and electron- 
ion collision strengths from Zhang & Pradhan (1995). We 
neglect fluorescent excitation by the UV continuum of the 
star. Indeed, as pointed out by Verner et al. (20 QBJ) in the 
case of the Orion Nebula, the 16 lowest levels of the Fe + 
ion remain collisionally dominated in the conditions un- 
der study (optically thin emission, n c = 1 — 10 8 cm -3 , 
T ~ 10 4 K). Thus, the line emissivity per Fe + ion only 
depends on electron density (n e ) and temperature (T e ). 

We have tested our simplified model by comparing our 
results with a more accurate 142-level computation, tak- 
ing into account the same physical processes (Pradhan, 
private communication). The comparison is shown in Fig. 
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Fig. 3. Diagnostic diagram of (n c ,T c ): Temperature sen- 
sitive ratio 0.8617 /tm/1.257 /jm versus electronic den- 
sity sensitive ratio 1.533 /um/ 1.644 /im. Dashed lines: 
n c varying from 1 to 10 6 cm -3 in factors of 10. Solid 
lines: T c varying from 5000 to 20 000 K (T c expressed in 
factors of 10 3 K). Error bars correspond to typical uncer- 
tainties of ± 10 % and ± 15 % in 1.533 /xm/ 1.644 /im and 
0.8617 jUm/1.257 /xm respectively. 



□ for the same four line ratios as in Pradhan & Zhang 
(1993), which probe all four lowest terms of Fe + . Our 
predicted line ratios are in very good agreement with 
those obtained from the 142-level model, for n c rang- 
ing from 10 to 10 8 cm -3 and T c ranging from 3000 to 
2 x 10 4 K. Differences do not exceed 1% except for transi- 
tions originating from the upper 4 P level (e.g. 0.8617 /im, 
0.9268 /on) and T c > 10 000 K, where our model under- 
estimates the 0.8617 /tm/1.257 /jm ratio by up to 5% for 
r±e < 10 5 cm -3 , and by — 20 % for n c > 10 5 cm -3 . These 
differences still remain below current observational uncer- 
tainties in line ratios. Our simplified 16-level model is thus 
sufficiently accurate for diagnostic purposes. 

In the following, we concentrate on the prominent 
optical and near-infrared lines. Most of the brightest 
transitions arise from the 4 D term, namely 1.257 /im, 
1.644 fim and 1.533 /mi (Fig.HJ. The decay 4 F 9/2 - 4 P 5/2 
(0.8617 /im) or e D 9 / 2 — 4 -F9/2 (5.339 /Ltm) lines also have 
significant emissivities (> one tenth of the 1.257 /xm line). 
The latter is inaccessible from ground-based telescopes 
but may be observable from space, with JWST for ex- 
ample. 

2.2. Derivation of (n e , T e ) 

Critical densities for the most intense [Fe Ii] lines, and 
for the main forbidden lines of [S 11] and [O 1] observed 
in jets, can be calculated with the generalized formula of 
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Table 1. Critical densities (n cr ) of upper levels for the 
prominent atomic transitions, calculated using Eq. (1) for 
T c = 5000 K and 10 000 K 



Osterbrock ((1989): 

n cr (i) = Sj<j Aij/Hjjti qij, (1) 

where Ay denotes Einstein spontaneous emission coeffi- 
cients and qij denotes collisional rate coefficients. The 
results are given in Table ^ an d show that the [Fe 11] 
lines investigated here have critical densities of a few 
10 4 cm" 3 (slightly higher than the red [S 11] AA6716.6731 
lines 1 ). They will therefore probe emitting regions with 
electronic densities of the same order. Line ratios in- 
volving transitions with similar excitation energies and 
similar ground states can be used to derive n c . For n c 
in the range 10 2 to 10 5 cm -3 , the best diagnostic is 
the [Fe 11] 1.644 /tm/1.533 /tm ratio. Indeed, this ratio 
involves two of the brightest near-IR [Fe Ii] lines and 
is little dependent on temperature and reddening (Fig. 
EJd). Note that the brightest two transitions at 1.644 and 
1.257 fiin originate from the same upper level. Their ratio 
(e(1.257 /tm)/e(1.644 /im) = 1.36) is therefore insensitive 
to the excitation conditions and can be used to derive the 
extinction Ay along the line of sight. 

Since most of the prominent near-infrared transitions 
have similar excitation energies, their ratio is not very sen- 
sitive to the gas temperature. Combining them with tran- 
sitions from the 4 P term, with higher excitation energy, 
could in principle lead to gas temperature determinations. 
In the near- infrared domain, the transitions at 1.811 /tm 
and 1.813 /im would be very good candidates. However, 
these lines lie in a region of poor atmospheric transmis- 
sion and in addition are blended with the stronger A F— A D 
1.810 /tm line. The line at 1.748 /tm is weaker and may 
be, at the present available spectral resolutions, strongly 
contaminated by the H2 1-0 S(7) transition (cf. Nisini et 
al. 2002). A few lines lying in the red part of the opti- 
cal spectrum can be used, the most prominent being at 
0.8617 /tm and 0.8892 /tm. In Fig. [31 we propose a di- 
agnostic diagram, combining the density sensitive [Fe 11] 
1.644 /tm/1.533 /im ratio with the temperature sensitive 
[Fe 11] 0.8617 /im/1.257 /im ratio, which provides a direct 
measure of (n e ,T e ). A similar diagram, including the red 
transition at 4 F3/2— 4 -P3/2 0.9471 /tm, is proposed in Nisini 

1 The critical density for the [S 11] A6731 line in Table [T] is 
markedly lower than the value ~ 10 4 cm -3 often quoted in the 
literature, derived from a two-level atom approximation. 
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Fig. 4. Diagnostic diagram of the Fe + /S + gas phase 
abundance. We plot the density sensitive ([Fe n] 
1.533 /um/1.664 /zm) ratio versus the depletion sensitive 
[Fe n] 1.257 /im/([S n] (1.029+1.032+1.034+1.037) pm 
ratio. Predicted line ratios are plotted as function of T e 
(solid lines) and n c (dashed lines) for Fc+/S + = 1.7 (corre- 
sponding to solar abundances, Savage & Sembach 11996(1 . 
Dotted and dash-dotted lines correspond to Fe + /S + = 
0.17 (1/10 of solar value). Observed line ratios for the sam- 
ple of HH objects studied by Nisini et al. <|2002fl . corrected 
for individual reddenings, are represented with their error 
bars (including Ay uncertainties). 



et al. H2()()2jl . With a flux ratio uncertainty of ± 10 % for 
the [Fe n] 1.644 /im/1.533 /im ratio, the diagnostic di- 
agram presented in Fig. [3] would allow determination of 
n ranging from 10 2 to 10 5 cm~ 3 with less than a fac- 
tor 4 uncertainty (total error). To achieve an uncertainty 
of ± 15 % on T c at T c < 10 000 K, an estimated error of 
± 15 % in the [Fe n] 0.8617 /mi/1.257 /jm ratio is required. 

2.3. Diagnostic of gas phase Fe abundance 

In the interstellar medium, condensation processes onto 
dust grains deplete the gas phase abundance of heavy el- 
ements such as Fe by a factor > 10 (Savage & Sembach 
1996). In jets from young stellar objects, sputtering pro- 
cesses in sufficiently fast shocks (see e.g. Jones et al. H994[l 
could destroy a significant fraction of the dust grains, re- 
plenishing the gas phase with heavy elements. Strong heat- 
ing at the jet base could also evaporate dust grains in the 
ejected gas. The gas phase abundance of Fe in HH flows is 
therefore an important parameter giving interesting indi- 
cations on the shock and thermal history of ejected matter. 
In addition, estimates of the Fe + gas phase abundance in 
HH flows are critically required to derive mass-loss rates 
from [Fe il] line fluxes. 



A first method to assess the iron gas phase depletion 
consists in comparing the Fe abundance to that of sul- 
fur, an element not significantly depleted onto dust grains, 
with a known gas-phase abundance of 1.86 x 10~ 5 (Savage 
& Sembach [1996). In young stellar jets, emission of these 
elements is dominated by [Fe n] and [S n] lines. However, 
the Fe + /S + ratio derived from these lines is probably very 
similar to the total Fe/S ratio. If the lines are formed 
in cooling flows behind shock waves, one expects Fe + /Fe 
< S+/S < 2Fe+/Fe (Beck-Winchatz et al. 1996), so that 
Fc+/S+ < Fc/S < 2Fc+/S+. If jet heating is instead dom- 
inated by ambipolar diffusion and UV radiation from the 
accretion shock, S and Fe should be both predominantly in 
their singly ionized form (see Garcia et al. 1200 la} , so that 
Fe/S = Fc+/S + . In either case, the Fe + /S + ratio provides 
a valuable lower limit to the actual Fe/S ratio, within a 
factor of 2 of the true value. 

This method was first applied to HH objects by Beck- 
Winchatz et al. {1994 1996), using a large sample of opti- 
cal [Fe n] lines. The Fe + /S + ratio is derived from the rel- 
ative intensities of [Fe n] and [S n] optical lines by means 
of a statistical equilibrium model, similar to that outlined 
in Sect. 2.1, given appropriate values of (n e ,T e ). The as- 
sumption of a uniform emission region is partly justified 
by the fact that the lines are excited under rather similar 
conditions (see Sect. 2.2). Inferred gas phase Fe + /S + ra- 
tios range from 1 to 4 when the [S n] AA6716,6731 lines are 
used for reference. When compared with the most recent 
estimate of the solar Fe/S ratio ~ 1.7, consistent with the 
local ISM (Savage & Semha,ch [TOSSI Sofia & Meyer HHJ), 
these values indicate no Fe gas phase depletion in most 
HH flows. Surprisingly, no correlation with the degree of 
excitation of the flow is observed. This analysis was how- 
ever limited to the case of optically revealed parts of the 
jets, i.e. usually far from the driving source. 

Recently, Nisini et al. I|2002|l proposed two other meth- 
ods to measure the Fe gas phase abundance in young stel- 
lar jets, relying only on near-infrared lines. The first one 
involves the [Fe n] 1.257 /zm/Pa/3 ratio, assuming that 
Pa/? is emitted under case B recombination. A drawback 
of this method is that, in addition to n G and T e , it requires 
knowledge of the hydrogen ionization fraction, x e , which 
is not available for all sources. Furthermore, the assump- 
tion that Pa/3 is formed in the same region as [Fe n], and 
that it is dominated by recombination, is generally not 
valid in shocks (cf. Bacciotti & Eisloffel 1999). A second 
method proposed by these authors is to compare the ob- 
served near-IR line ratios [Fe n] 1.257 /zm/[C i] 0.98 /jm 
and [Fe n] 1.257 ^m/Pa/3 with predictions from models 
of fast dissociative shocks. In the sample of objects stud- 
ied by Nisini et al. |2002), the inferred gas phase Fe/H 
ranges from 30 % to 70 % of the solar photospheric value 
((Fe/H) = 3.2 x 10" 5 , Savage & Sembach ESS- These 
results depend on the assumption that the carbon abun- 
dance in the shock remains constant at 2.3 x 10~ 4 , i.e. 64 % 
of the solar value (Hollenbach & McKee 119891 Savage & 
Sembach I1996[) . while Fe is returned to the gas phase. A 
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solar gas-phase carbon abundance would increase their es- 
timated Fe/H by typically a factor 1.5. 

We propose here a third approach to estimate the 
gas phase Fe abundance of HH flows from near-IR ob- 
servations, independent of both x e and shock models. 
The method is simply an extension to the near-IR do- 
main of that introduced by Beck-Winchatz et al. {1994) 
in the optical range. As described above, it will yield 
a lower limit to the Fe/S ratio. We propose to use the 
ratio of the [Fe n] 1.257 (im line to the [S n] 1.029 + 
1.032 + 1.034 + 1.037 /im lines (hereafter referred to as 
[S n] 1.03 /im). We predict comparable integrated flux 
in these lines in the conditions under study and, indeed, 
these near-IR [S n] lines have been detected by Nisini et al. 
(2002) in their sample of HH objects. We have chosen the 
[Fe n] 1.257 /im line because it is one of the brightest tran- 
sitions and it minimizes reddening uncertainties relative to 
[S n] 1.03 /im. To compute the emissivity per ion for S + , 
we construct a model including the 5 low-lying metastable 
levels. Transition probabilities are taken from Fritzsche et 
al. I|1999|) . collision strengths from Cai & Pradhan Q1993[l . 
We also use the [Fe n] 1.533 /im/1.644 /im ratio, as an 
indicator of n . Fig. 0] presents the resulting diagnostic 
diagram, which plots [S n] 1.03 /im/[Fe n] 1.257 /im ver- 
sus [Fe n] 1.533 /im/1.644 /im, and yields the gas phase 
Fe + /S + ratio provided T c is known. It can be seen that 
near-IR [Fe n] / [S n] line ratios are more sensitive to tem- 
perature than the optical ratios (see Table 2 of Beck- 
Winchatz et al. 11994), though they are less sensitive to 
reddening. 

We plot in Fig. 0] the observational data from Nisini 
et al. (J2D02). Temperatures are not available for these 
sources, but a range of 7000 to 15 000 K can be safely 
assumed, based on the predicted temperature of the low 
excitation forbidden line emitting region in planar shock 
models on the one hand (Hartigan et al. I1994fl . and max- 
imum derived temperatures in optically visible HH flows, 
on the other hand (Bohm & Solf 1990). We can derive 
an absolute lower limit to the Fe + /S + gas phase ratio, 
and hence to the Fe/S ratio, by assuming T c = 7000 K. 
These limits are: 0.25 ± 0.1 for HH 111H, 0.39 ± 0.1 for 
HH 241A, 0.55±0.2 for HH 111F, 0.83±0.15 for HH 120H, 
and 1.2 ± 0.25 for HH 240A, and range from ~ 15 % to 
70 % of the solar value (Fe/S = 1.7). These values are 
remarkably similar to the results found by Nisini et al. 
( 2002J) from shock models, with the same ordering for the 
5 individual sources. However, we stress that the above 
estimates are lower limits. With T e ~ 10 4 K, the data 
are compatible with no depletion in most sources, except 
HH 241 A and HH 111H (50 % depletion). In the follow- 
ing, we will thus make the assumption that all iron is in 
the gas phase. However, we point out that to derive iron 
depletion and jet mass-loss rates to better than a factor 
~ 5 from [Fe n] line fluxes, T estimates will be critically 
required. 



3. Predictions for cold MHD disk winds 

In this section, we use our simplified Fe + model to pre- 
dict near-IR [Fe n] emissivities for the class of self-similar 
MHD cold disk wind models developed by Ferreira 0997 ) , 
for which fully consistent thermal solutions have been re- 
cently computed (Garcia et al. I2001afl . We extend here 
the work presented in Garcia et al. I|2001bj) . which was 
restricted to predictions in the optical domain. We first 
briefly describe the model and compare predicted emissiv- 
ities in the near-IR [Fe n] lines with those of the prominent 
optical [O i] A6300 and [S il] A6731 emission lines. We then 
compare these predictions with recent observations in the 
near-IR domain of the stellar jets from L1551TRS5 and 
DG Tau. 



3.1. Wind model 

The MHD accretion-ejection structures of Ferreira (|l997) 
represent an improvement on the centrifugal MHD disk 
winds first calculated by Blandford & Payne l|1982fl . Both 
assume that (1) matter is ejected along a large scale bipo- 
lar magnetic field threading through the accretion disk, 
(2) jet enthalpy is negligible for accelerating the flow 
(the wind is "cold"), and (3) the structure is steady- 
state, axisymmetric, and self-similar to the disk radius. 
However, the solutions of Ferreira Ijl997|l also describe 
self-consistently the accretion flow in the underlying resis- 
tive Keplerian disk, including all relevant dynamical terms 
that govern the wind mass-loading at the slow-point (see 
Ferreira & Pelletier ll995|) . This global description imposes 
additional constraints on the self-similar structure, which 
is then specified by only three dimensionless parameters: 
(1) £ = dlogM acc /dlogT, the ejection efficiency parame- 
ter, which is related to the total ejection/accretion ratio 
by M e j ec /Macc — £ x ln(r e /n), where r e and r, are the 
outer and inner radii of the MHD disk- wind structure, (2) 
e = Ii/tq , where h is the disk scale height at the cylindrical 
radius ro, and (3) a m = v m /VAh, where v m is the required 
turbulent magnetic diffusivity and Va is the Alfvn velocity 
on the disk midplane. 

Following Cabrit et al. (19.93); we wm adopt e = 0.1 (as 
estimated in HH 30 by Burrows et al. [19 9 51 and a m = 1. 
Solutions that extend far from the Alfvn surface are then 
found for £ between 0.005 and 0.012 (Ferreira 1997). We 
concentrate here on the solution with £ = 0.01 (Model A 
in Cabrit et al. I1999f) . which is found to best reproduce 
the collimation properties of T Tauri microjets (Cabrit et 
al. 119991 Dougados et al. I2000j). As a result of toroidal 
magnetic fields, the streamlines slowly recollimate beyond 
z/ro ~ 1500 (where ro is the initial radius of the streamline 
in the disk), and terminate at z/tq ~ 8000, due to strong 
rcfocusing towards the axis. One expects that a more grad- 
ual recollimation would occur in reality, due to additional 
pressure on the flow axis from open stellar magnetic field 
lines. 



6 



N. Pesenti et al.: Near-IR [Fe n] emission. 



3.2. Predicted [Fe n] emissivities 

The thermal and ionization structure of the MHD disk 
wind solutions of Ferreira H1997[) was computed a poste- 
riori by Garcia et al. f 2 001a|l for an atomic gas of solar 
composition, along 11 streamlines anchored in the disk 
at radii ranging from ro = 0.07 AU, the typical disk 
corotation radius for a T Tauri star, to 1 AU, where 
molecules should start to form. The dominant heating 
mechanism is ambipolar diffusion, while adiabatic expan- 
sion dominates the cooling. We will use their results for 
a central stellar mass M* = 0.5 M Q and accretion rates 
Mace = 10~ 5 — 10~ 6 M yr _1 . We compute the emissivi- 
ties in the [Fe n] near-IR lines along each streamline with 
the 16-level atom model described in Section 2.1, using 
the n c , T and Fe + /Fe ratios provided by the thermal so- 
lution. The emissivity grid is projected onto the plane of 
the sky to construct emission maps and long slit spectra. 
Maps are then convolved by a two-dimensional Gaussian 
beam, representative of the instrumental spatial and spec- 
tral resolutions. 

We first compare predicted near-IR [Fe n] emissivities 
along the flow with those predicted for the main opti- 
cal emission lines of [S n] and [O i]. We plot in Fig. [S] 
the variation of emissivity with altitude along the inner- 
most and outermost streamlines for the [Fe n] 1.644 jitm, 
[S II] A6731 and [O i] A6300 lines. Both streamlines follow 
the same behavior in all 3 lines: emissivity first increases 
due to the central rise in temperature, the plasma then 
reaches a temperature plateau (~ 10 4 K) and the drop 
in n e , due to the widening of the jet, produces the subse- 
quent emissivity decrease. The final increase in the inner 
streamline at Z > 200 AU is due to compression in the 
recollimation zone. Since solutions have a self-similar ge- 
ometry, maximum emissivity is reached at a lower altitude 
(Z < 1 AU) in the inner streamline as compared to the 
outer one (Z - 10 AU). 

There are however important differences in emissiv- 
ity profile between the 3 lines close to the central source. 
The upper level energies of near-IR [Fe n] lines (E/A;~ 
10 000 K) are lower than for [S n] A6731 and [O i] A6300 
lines (E/k - 20 000 K), so that the [Fe n]/[S n] and 
[Fe n]/[0 i] ratios increase dramatically at the moderate 
electron temperatures < 7000 K prevailing at the wind 
base. [Fe n] lines therefore appear to be a choice tracer for 
the cooler, dense regions at the wind base, typically below 
a few AUs. In addition, these near-IR lines would allow 
to better trace extincted regions at the base of embedded 
young stellar jets. 

The AMBER/VLTI instrument, with a spatial resolu- 
tion of ~ 1 mas (0.14 AU in Taurus) in the near-IR, will 
be ideally suited for such studies of innermost jet regions. 
We plot in Fig. the predicted AMBER/VLTI map in 
[Fe n] 1.644 fim for our fiducial cold disk wind model, with 
M acc = 10~ 6 M Q yr _1 and inclination angle to the line of 
sight of 45°. We also show as white contours the ratio of 
the [Fe n] 1.644 ftm to the [S il] A6731 surface brightnesses. 
As expected, the [Fe n] line is several times stronger than 
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Fig. 5. Top: Variation of n c and n# (bold curve) along 
the inner (ro = 0.07 AU, full lines) and outer (ro = 1 AU, 
dashed lines) streamlines of the fiducial cold disk wind 
solution (£ = 0.01, Macc = 10~ 6 Moyr" 1 , M* = 0.5 M Q ) 
as a function of the altitude Z above the disk. Middle and 
Bottom: Variation of line emissivities along the inner and 
outer streamlines for the 3 emission lines: [Fe n] 1.644 /im, 
[S n] A6731 and [O i] A6300. The variation of T e is also 
shown in bold curve. 



the [S n] one on scales < 2 AU above the disk, due to 
the drop in T e . At this spatial resolution, the [Fe n] map 
is dominated by the innermost streamline, which is fully 
resolved, and shows a clear central hole. Hence, such ob- 
servations could provide direct constraints on the shape 
(Z(r)) of the streamlines in the MHD jet solution and on 
the innermost radius at which ejection occurs. 



3.3. Comparison with L1551-IRS5 and DG Tau 

In this section, we discuss constraints on the ejection 
mechanism brought by recent observations in near-IR 
[Fe II] lines of the L1551-IRS5 and DG Tau jets. We com- 
pare these observations with predictions from our fiducial 
cold disk wind model and with constraints derived from 
previous studies of T Tauri jets in the optical domain. 



N. Pesenti et al.: Near-IR [Fe n] emission. 



7 




Fig. 6. Predicted [Fe n] 1.644 /im emission map of 
the fiducial cold disk wind solution with £ = 0.01, 
M aC c=10~ 6 Moyr -1 and inclination to the line of sight 
i of 45°. The field of view (0.04") and the Gaussian beam 
(1 mas, i.e. 0.14 AU) correspond to the AMBER/VLTI 
instrument characteristics. Filled contours start at the 
maximum surface brightness (2.3 10 -4 Wm _2 st _1 ) and 
decrease by factors of 2. White contours correspond to 
[Fe n] 1.644 /mi/[S n] A6731 ratios of 10 and 1 (from left 
to right). The shape of the innermost streamline, assumed 
here anchored at 0.07 AU, is fully resolved and would set 
strong constraints on proposed MHD models. 



3.3.1. Morphology and kinematics 

One of the best-known young stellar jets is associated with 
the embedded binary source L1551-IRS5. It was first de- 
tected in the optical by Cohen et al. (1982). High reso- 
lution images show two extended ridges, identified as two 
separate jets with differing line of sight velocities (Fridhmd 
& Liseau 115551 Itoh et al. l2UUU|l . Pyo et al. (|2UU2jl present 
velocity-resolved [Fe n] 1.644 /im spectra towards the base 
of the northern jet. These observations provide detailed 
morphological and kinematical information on the inner 
500 AU of the jet, which can be directly compared to pre- 
dictions from the cold disk wind model investigated here. 

We adopt for the model an accretion rate of M acc = 
10 -5 Mq yr _1 close to the estimated value for the L1551- 
IRS5 source of 6 x 10" 6 Moyr" 1 (Momose et al. EHHJ), 
an inclination angle to the line of sight of 45° (Pyo et al. 
20021, an d a central stellar mass of 0.5 M Q , i.e. roughly 
half the total mass of the binary system ~ 1.2 M Q es- 
timated from orbital motions (Rodriguez et al. [2003 ) . 
Predicted emission maps and position-velocity (PV) dia- 
grams are convolved with Gaussian beams of widths 0.3" 
and 60 kms -1 , corresponding to the instrumental spatial 
and spectral resolutions of the L1551 observations. 



Fig. Ek, b,c present a comparison of observed and pre- 
dicted maps in the [Fe n] 1.644 /im line. We also show 
in Fig. 01 synthetic maps taking into account the strong 
extinction observed towards L1551-IRS5. We artificially 
dimmed the predicted maps, before convolution with the 
instrumental beam, by a uniform screen of A v = 20 mag 
extending out to 1", following the extinction measure- 
ments of Itoh et al. PTOHjl . 

The synthetic images show that our cold disk wind 
model reproduces quite well the general observed mor- 
phology (although a higher extinction towards the source 
would be needed to hide completely the base of the jet). 
Interestingly, the sharp decrease in extinction at 1" pro- 
duces an emission peak close to the observed peak PHK 1. 
We also note that refocussing onto the jet axis in our 
model produces a peak at a distance comparable to the 
distance of knot PHK 2. In synthetic PV diagrams, the 
refocusing of streamlines produces a strong emission ridge 
beyond 2" from the star, as the electronic temperature and 
density are locally enhanced by compression. The appar- 
ent deceleration in this ridge is due to the self-similarity 
of the solution: Streamlines of lower speeds are launched 
from larger disk radii, and thus recollimate farther from 
the star. A dashed curve in Fig. 0: shows the locus of the 
points where streamlines start to recollimate. Predictions 
downstream from this curve are only indicative, since gas 
pressure effects (neglected in the MHD solution) should 
eventually counteract compression, and start to have an 
important effect on the dynamics. 

Fig. plots the synthetic line profile in the region 
before jet recollimation occurs. It shows that our cold 
disk wind model reproduces quite remarkably the main 
kinematical properties of the central regions of the L1551 
jet: the observed velocity full width at zero intensity of 
300 kms -1 (350 kms -1 in the model) as well as the exis- 
tence of two velocity components: a high velocity compo- 
nent (HVC) and an intermediate- velocity component (IVC 
2 ) with centroid velocities —300 kms -1 (—360 kms -1 in 
the model), and —100 kms -1 (—160 kms -1 in the model) 
respectively. We stress in particular that a single cold disk 
wind model can fully account for the global kinematical 
properties, in particular the velocity extent of the pro- 
file, due to the range of streamlines involved. Streamlines 
originating from the disk inner truncation radius rapidly 
collimate, producing a fast narrow jet, while outer stream- 
lines produce a wider angle slower wind accounting for the 
low-velocity part of the profile. As already pointed out 
in Cabrit et al. (|1999|) and Garcia et al. (|2001bjl . a two 
component wind model, as earlier suggested by Kwan & 
Tademaru (^)95), may not be necessary to account for 
both observed velocity peaks (HVC and IVC). 

We notice however one important discrepancy between 
model and observations, apart from the apparent deceler- 

2 We use this denomination to differentiate the IVC from the 
true low velocity component observed in classical T Tauri stars 
which have much lower centroid velocities ( — 10 kms -1 ) and 
are spatially unresolved. 
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Fig. 7. [Fc n] 1.644 /im emission maps (Left) and 
position- velocity diagrams along the jet (Right). (a),(b): 
Observations of the L1551 jet from (Pyo et al. 2002) with 
a seeing size of 0.3" and a velocity resolution of 59 kms . 

(c) : synthetic maps predicted for a cold disk wind model 
with I = 0.01, M acc = 10~ 5 M Q yr-\ i = 45°, convolved 
with the beam sizes of the L1551 observations. The cross 
and dot-dashed lines locate the position of IRS5 VLA as 
in the top panel. The dashed line traces the locus where 
streamlines start to recollimate (beyond which our predic- 
tions are only approximative; see text for more details). 

(d) : Synthetic maps including the effect of an extinction 
screen towards L1551-IRS5. (e): Line profile extracted 
from the predicted PV diagram (Panel (c)) at the posi- 
tion of PHK1. 



ation of the high-velocity emission ridge discussed earlier 
and due to the refocusing of streamlines towards the axis 
in the model. The relative intensity of the IVC and HVC 
at large distances is not well reproduced by the cold disk 
wind model. In the L1551 observations, the intermediate- 
velocity component dominates the emission out to 3" 
from the star, while in the model the IVC emission drops 
rapidly beyond 0.5". We note the possibility that the 
position-velocity diagram of Pyo et al. ( 2002) is contami- 
nated by the southern L1551 jet (with line of sight veloc- 
ities of ~ —60 kms -1 , Hartigan et al. I2000f) . especially at 
intermediate- velocity towards the PHK 1 knot where the 
two jets appear to merge in the images. However, follow- 
ing Pyo et al. ( 2 0712^1 . we estimate this possibility unlikely 
at distances > 1" where the southern jet is several times 
fainter than the northern one. 

We now apply the same type of comparison to DG Tau, 
an "extreme" T Tauri star with an infrared excess inter- 
mediate between Class I and Class II, and a particularly 
bright and well-studied optical jet (Lavalley-Fouquet et al. 
2000] Dougados et al. 123001 Bacciotti et al. l2T)tI))l . A direct 
comparison with models is possible since the jet inclina- 
tion is known from proper motions, and high-resolution 
velocity-resolved [Fe n] 1.644 /im spectra were recently 
published by Pyo et al. H2()()3|l . In DG Tau, ejection signa- 
tures have been shown to be highly time variable close to 
the central source: Large changes in the measured velocity 
distribution, covering the range from to —400 kms -1 , 
are observed at distances < 0.5-1" on timescales of a few 
years (Solf H99Tjl . However, beyond these distances, ter- 
minal velocities varied by less than 20 % over 10 years: 
-200 kms -1 at d = 3-3.5" observed in Solf & Bohm 
ipg5|) . -210 to -190 kms -1 at d = 2.5-4" in Lavalley- 
Fouquet et al. PTOjl . -240 kms" 1 at d > 1.2" in Pyo et 
al. ( 2003). Internal shocks due to variability of the ejec- 
tion velocity are expected to dissipate quickly along the jet 
(Raga & Kofman 1992), leaving the average flow mostly 
unperturbed at large distances. Indeed, derived shock ve- 
locities in the knots of the DG Tau jet do not exceed 50- 
70 kms -1 (Lavalley-Fouquet et al. 2000). A direct com- 
parison of the properties of the jet on scales > 100 AU 
with the predictions of stationary models appears there- 
fore justified. 

In Fig. |SJ we compare the long-slit spectrum from Pyo 
et al. <|2003|) with our model predictions, convolved by the 
appropriate resolution (0.16" and 30 kms -1 ). We took an 
accretion rate of 10~ 6 M Q yr _1 , close to the estimated 
value of 2 x 10~ 6 M Q yr" 1 through the disk (Hartigan 
et al. I1995|l and an inclination angle of 45°, as estimated 
from the proper motions (Dougados et al. 2000. Pyo et al. 
2003). The model again explains well the presence of two 
velocity components, and the general centroid velocities: 

the IVC peaks at 60 to -130 kms -1 in DG Tau (-100 

to -140 kms -1 in the model), the HVC peaks at -200 
to -240 kms -1 in DG Tau (-360 kms" 1 in the model). 
However, we note that the agreement for the HVC cen- 
troid is less good than in L1551-IRS5, in the sense that 
the model now overestimates the HVC velocity by a fac- 
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Fig. 8. [Fe n] 1.644 /mi position-velocity diagrams. 
Contours are drawn with equal intervals in logarithmic 
scale. Top: Observations of the DG Tau jet (Pyo et al. 
2003J with an angular resolution of 0.16" and a veloc- 
ity resolution of 30 kms" 1 . Bottom: Predicted diagram 
from the cold disk wind model A (Garcia et al. [2001b) 
with £ = 0.01, M acc = 1Q- 6 M Q yr-\ % = 45°, convolved 
with the beam sizes of DG Tau observations. 

tor 1.5. Another discrepancy, also encountered in the case 
of L1551-IRS5, is that the IVC rapidly becomes weaker 
than the HVC in our model (at Z ~ 0.2"), while in the 
observations of Pyo et al. H2003[) . the IVC dominates over 
the HVC much further out from the star (until Z ~ 0.6"). 

3.3.2. Density 

Itoh et al. 1|2000[) derive the variation of extinction and 
n e along the L1551 jet from JH-band spectroscopic ob- 
servations, using the [Fe ii] 1.257 /im/1.644 /im and 
1.644 /im/1.533 /im line ratios (see section 12.2(1 . The 
derived electron densities decrease from 10 6 cm~ 3 (at 
d=2" = 280 AU from IRS5) to 10 4 cm- 3 (d = 4" 



from IRS5). For our fiducial cold disk wind model (with 
M acc = 10~ 5 M Q yr -1 ), we compute n c from [Fe n] 
1.644 /im/1.533 /im line ratios integrated across the jet. 
On the same spatial scales the predicted value is ~ 
10 3 cm -3 , i.e. 10 to 1000 times smaller than observed. 
We note that even assuming fully ionized gas (i.e. extra 
heating mechanisms: shocks or turbulent dissipation), the 
density in our model would still remain too low: all stream- 
lines have ti h < 10 6 cm~ 3 beyond Z = 100 AU for M acc = 
10~ 5 M yr" 1 (see Fig.EJwith M acc = 10~ 6 M yr" 1 , and 
note that n# oc M acc ). This density deficiency has already 
been observed in the central regions of T Tauri micro- 
jets (Lavalley-Fouquet et al. 120001 Garcia et al. I2001bjl . 
The near-IR observations of the L1551 jet show that self- 
similar cold disk winds also fail to reproduce the observed 
electronic densities in younger embedded flows on spatial 
scales of a few 100 AUs. 

3.3.3. Summary 

While the near-IR [Fe n] profiles predicted by the cold disk 
wind model reproduce rather well the two velocity compo- 
nents (HVC and IVC) observed in stellar jets, they appear 
to suffer from the same deficiency of intermediate- velocity 
emission (as compared to high-velocity) already noted in 
optical studies of the inner jet regions in T Tauri stars 
(Cabrit et al. IT533I Garcia et al. I2001bjl . More efficient 
heating at the base of the jet appears needed to enhance 
the brightness of the IVC component compared to the 
HVC. In addition, the predicted terminal velocity for the 
HVC appears in very good agreement (within ~ 20 %) 
with the observations in the Class I source L1551-IRS5, 
but exceeds the observed speed of the HVC by a factor 
~ 1.5 in the Class I/II source DG Tau. We note that the 
asymptotic poloidal velocity predicted by the disk wind 
models can be expressed by: v po ;,oo — «k,ov / 2A — 3, where 
Vk,o is the Keplerian speed in the disk midplane at the 
cylindrical launching radius ro and A ~ r^Ao the mag- 
netic lever arm (ta is the radius of the magnetic surface 
at the Alfven point). Lower terminal poloidal velocities 
can be achieved by decreasing either vk,q (i.e. increasing 
ro) or A. Larger launching radii would lead however to 
even lower total jet densities, increasing the discrepancy 
with the observations. Therefore, smaller magnetic lever 
arms (typically A < 25 as compared to our disk wind solu- 
tion where A ~ 50) are required to reproduce the observed 
terminal speed in DG Tau. 

4. Conclusions 

Thanks to a simplified 16-level atomic model, we investi- 
gated the potential of near-IR [Fe n] lines as diagnostics 
of physical conditions in the inner regions of stellar jets. 
We presented diagnostic diagrams that provide determi- 
nations of n c in the range 10 2 to 10 5 cm" 3 , T c below 
20 000 K, and a lower limit to the Fe/S ratio (based on a 
near-infrared [Fe n]/[S n] line ratio). Using observational 
data from Nisini et al. (2002) in several HH flows, and as- 
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suming T c ~7000 K, we deduce lower limits to the iron gas 
phase abundance comparable to those found from molec- 
ular shock models by Nisini et al. ( 2002}, i.e. on average 
30% of the solar value. For T c ~ 10 4 K, the derived Fe 
abundances would be roughly solar in most HH objects, 
except HH 111H and HH 241 A. Hence, we confirm earlier 
work indicating an efficient dust destruction mechanism 
in stellar jets. However, we cannot yet conclude on the 
nature of the grain destruction process, i.e. shocks versus 
evaporation. Calculations of grain destruction in hydrody- 
namic shocks have concentrated so far on the low-density 
interstellar medium (Jones et al. H994[l and similar compu- 
tations with preshock densities ~ 10 3 — 10 4 cm~ 3 typical 
of HH flows are needed. Such calculations currently exist 
only in the case of low-ionization magneto-hydrodynamics 
shocks with ion-neutral decoupling, but this type of shocks 
may not occur at the ionization level of a few % observed 
in optical jets. High angular resolution observations of the 
jet base in [Fe n] and [S n] lines, e.g. with NAOS/CONICA 
on the VLT, are also needed to determine whether the 
grains are destroyed in the vinicity of the exciting source 
or further out. 

In the second part of this article, we extended the work 
presented in Garcia et al. i|2001bj> to predictions in the 
near-infrared domain, by applying our Fe + emissivity cal- 
culations to a cold disk wind model with a self-consistent 
thermal and ionization structure (Ferreira 119971 Garcia 
et al. I2001a|l . On large scale (a few 100 AUs), predicted 
emission maps are very similar in near-IR [Fe II] lines and 
optical [S n] lines. However, on scales of 1 AU soon ac- 
cessible with the AMBER/VLTI instrument, [Fe il] lines 
appear more efficient to trace colder and more extincted 
regions very close to the driving source. In particular, 
AMBER/VLTI observations in [Fe n] lines are predicted 
to resolve the innermost streamline and inner hole ex- 
pected for a cold disk wind. 

Comparing synthetic maps with recent observations of 
the L1551-IRS5 and DG Tau jets in near-infrared [Fe n] 
lines, we find that our model is in good agreement with the 
global kinematics, and notably the presence of two veloc- 
ity components: a intermediate- velocity component (IVC) 
at —100 kms -1 (—100 to —150 kins" 1 in the model) and a 
high-velocity component (HVC) at —230 to —300 kms -1 
(—360 kms -1 in the model). However, the predicted ve- 
locity for the HVC tends to be too high (by 20 % only in 
L1551 IRS5, but by a factor 1.5 in DG Tau). Furthermore, 
the model predicts insufficient emission at intermediate- 
velocity, compared to the HVC, and also insufficient gas 
density in L1551. Similar failures were encountered when 
comparing the cold disk wind model with observations of 
jets from T Tauri stars in the optical domain (Garcia et al. 
2001b). MHD disk wind models with entropy injection at 
the base of the jet (Casse & Ferreira PiOOOf) appear promis- 
ing to solve these problems as they could provide denser 
solutions, with a slower HVC and more emission from the 
IVC (originating from the warm wind base). Calculations 
of observational predictions for this class of "warm disk 



winds" are in progress and will be the subject of a future 
paper. 

Acknowledgements. We are very grateful to Anil Pradhan for 
communicating unpublished results from his 142-level code and 
for useful suggestions. We would like to thank the referee, K.- 
H. Bohm, for his helpful comments. 

References 

Bacciotti, F., & Eisloffel, J. 1999, A&A, 342, 717 

Bacciotti, F., Mundt, R., Ray, T. P., et al. 2000, ApJ, 537, L49 

Beck-Winchatz, B., Bohm, K. H., & Noriega-Crespo, A. 1994, 

PASP, 106, 1271 
Beck-Winchatz, B., Bohm, K. H., & Noriega-Crespo, A. 1996, 

AJ, 111, 346 

Blandford, R. D., & Payne, D. G. 1982, MNRAS, 199, 883 
Bohm, K.-H., & Solf, J. 1990, ApJ, 348, 297 
Burrows, C. J., Stapelfeldt, K. R., Watson, A. M., et al. 1996, 
ApJ, 473, 437 

Cabrit S., Edwards, S., Strom, S. E., & Strom, K. M. 1990, 
ApJ, 354, 687 

Cabrit S., Ferreira, J., & Raga, A. C. 1999, A&A, 343, L61 
Cabrit S., 2002 in Star Formation and the Physics of Young 

Stars, eds. J. Bouvier and J. P. Zahn, EAS Publications 

Series, Vol. 3, p. 147 
Cai, W, & Pradhan, A. K. 1993, ApJS, 88, 329 
Casse, F., & Ferreira, J. 2000, A&A, 361, 1178 
Cohen, M., Bieging, J. H., & Schwartz, P. R. 1982, ApJ, 253, 

707 

Dougados, C, Cabrit, S., Lavalley, C. & Mnard, F. 2000, A&A, 
357, L61 

Ferreira, J., & Pelletier, G. 1995, A&A, 295, 807 

Ferreira, J. 1997, A&A, 319, 340 

Fridlund, C. V. M., & Liseau, R. 1998, ApJ, 499, L75 

Fritzsche, S., Fricke, B., Geschke, D., Heitmann, A., & 

Sienkiewicz, J. E. 1999, ApJ, 518, 994 
Garcia, P. J. V., Ferreira, J., Cabrit, S., & Binette, L. 2001, 

A&A, 377, 589, Paper I 
Garcia, P. J. V., Cabrit, S., Ferreira, J., & Binette, L. 2001, 

A&A, 377, 609, Paper II 
Hartigan, P., Morse, J., & Raymond, J. 1994, ApJ, 436, 125 
Hartigan, P., Edwards, S., & Ghandour, L. 1995, ApJ, 452, 736 
Hartigan, P., Morse, J., Palunas, P., Bally, J., & Devine, D. 

2000, AJ, 119, 1872 
Hirth, G. A., Mundt, R., & Solf, J. 1997, A&AS, 126, 437 
Hollenbach, D. & McKee, C. F. 1989, ApJ, 342, 306 
Itoh, Y., Kaifu, N., Hayashi, M., et al. 2000, PASJ, 52, 81 
Jones, A. P., Tielens, A. G. G. M., Hollenbach, D. J., & McKee, 

C. F. 1994, ApJ, 433, 797 

Kwan J., & Tademaru, E. 1995, ApJ, 454, 382 
Lavalley-Fouquet, C, Cabrit, S., & Dougados, C. 2000, A&A, 
356, L41 

Momose, M., Ohashi, N., Kawabe, R., Nakano, T., & Hayashi, 

M. 1998, ApJ, 504, 314 
Nisini, B., Caratti o Garatti, A., Giannini, T., & Lorenzetti, 

D. 2002, A&A, 393, 1035 

Nussbaumer, H., & Storey, P. J. 1988, A&A, 193, 327 
Osterbrock, D. E. 1989, Astrophysics of Gaseous Nebulae and 

Active Galactic Nuclei, Univ. Science Books, Mill Valley 
Pradhan, A. K., & Zhang, H. L. 1993, ApJ, 409, L77 
Pyo, T., Hayashi, M., Kobayashi, N., et al. 2002, ApJ, 570, 724, 

Erratum, http : //www. subarutelescope . org/ staff /pyo/ 

L1551paper/Erratum . pdf 



N. Pesenti et al.: Near-IR [Fe n] emission.. 



11 



Pyo, T. -S., Kobayashi, N., Hayashi, M., et al. 2003, ApJ, 590, 
340 

Raga, A. C, & Kofman, L. 1992, ApJ, 386, 222 

Reipurth, B., Yu, K. C, Heathcote, S., Bally, J., Rodriguez, 

L. F. 2000, AJ, 120, 1449 
Rodriguez, L. F., Curiel, S., Canto, J., et al. 2003, ApJ, 583, 

330 

Savage, B. D., & Sembach, K. R. 1996, ARA&A, 34, 279 
Shang, H., Shu, F. H., & Glassgold, A. E. 1998, ApJ, 493, L91 
Shang, H., Glassgold, A. E., Shu, F. H., & Lizano, S. 2002, 

ApJ, 564, 853 
Sofia, U. J., & Meyer, D. M. 2001, ApJ, 554, L221 
Solf, J., & Bohm, K. H. 1993, ApJ, 410, L31 
Solf, J. 1997, in Herbig-Haro Flows and the Birth of Stars, IAU 

Symposium 182, Eds.: B. Reipurth & C. Bertout, Kluwer 

Academic Publishers, p. 63-72. 
Verner, E. M., Verner, D. A., Baldwin, J. A., Ferland, G. J., 

& Martin, P. G. 2000, ApJ, 543, 831 
Zhang, H. L., & Pradhan, A. K. 1995, A&A, 293, 953 



